Tunable coupling scheme for flux qubits at the optimal point 
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We discuss a practical design for tunably coupling a pair of flux qubits via the quantum inductance 
of a third high-frequency qubit. The design is particularly well suited for realizing a recently 
proposed microwave-induced parametric coupling scheme. This is attractive because the qubits can 
always remain at their optimal points. Furthermore, we will show that the resulting coupling also 
has an optimal point where it is insensitive to low-frequency flux noise. This is an important feature 
for the coherence of coupled qubits. The presented scheme is an experimentally realistic way of 
carrying out two-qubit gates and should be easily extended to multiqubit systems. 
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I. INTRODUCTION 

Superconducting qubits 1 " have received a lot of atten- 
tion during the last few years. The observed coherence 
times have been improved from the initial ns 2 range to 
the ^gMiSiSii range or so with the invention of so-called 
optimal bias points. That is, it is possible to set the 
biases of a qubit such that the energy difference of the 
utilized computational states do not depend on bias value 
to, say, first order. This makes the system at hand insen- 
sitive to harmful low-frequency noise. It seems that if an 
optimal bias point exists then the qubits should always 
be biased there 5 ^. 

However, isolated qubits are not usable for quantum 
computing. Some type of controllable coupling between 
individual qubits is desirable. Typically tunable cou- 
pling sc hemesi2ii2iiiii2iiiii2ii£il& require moving away 
from the optimal bias point. This is because the "natu- 
ral" coupling — inductive for flux qubits and capacitivc 
for charge qubits — is off-diagonal (a x ® a x ) at the opti- 
mal point and affects the qubit dynamics only in second 
order provided that the qubits are detuned. Therefore, 
constant moderate coupling can often be neglected at the 
optimal point. Recently, however, two promising schemes 
for coupling superconducting qubits (flux or charge) at 
the optimal point have emerged. In the so-called FLIC- 
FORQ coupling scheme^ two a x <g> er^-coupled detuned 
qubits can be made to interact by applying resonant mi- 
crowaves on the qubits at a power such that the sum 
of the Rabi frequencies of the individual qubits matches 
the detuning between the qubits. This causes a kind of 
a second resonance in the interaction picture, which can 
be used to realize a universal two-qubit gate. However, 
the scheme is challenging to realize currently with flux 
qubits since achieving small enough detunings in the fab- 
rication is hard and they easily exceed achievable Rabi 
frequencies. Nevertheless, because of the above reasons, 
Bertet et. al. 1 - suggested an alternative parametric cou- 
pling scenario which relies on the modulation of the cou- 



pling energy itself at the sum or the difference frequency 
of two detuned qubits. Liu et. al. 16 suggested earlier a 
related microwave coupling scheme which does not unfor- 
tunately work at the optimal point. The gained benefit 
in the scheme of Bertet et. al. compared to FLICFORQ 
is that the method is in principle applicable to qubits 
with a larger detuning and it may be possible to carry 
out gates faster. This scheme of course has an analogy 
for charge qubits too, but was suggested primarily for 
flux qubits. We also concentrate on flux qubits in this 
paper. The problem then is what kind of a physical cou- 
pling element to use to achieve the parametric control 
of the coupling energy. In the original paper— a method 
utilizing a current biased DC-SQUID— was considered to 
exemplify the scheme. This requires the modulation of 
the bias current at a high frequency. Moreover, the DC 
current bias of the SQUID must be near switching cur- 
rent and the flux near half flux quantum. This condition 
however contradicts with the optimal current bias con- 
dition of individual qubits 5 since the bias current noise 
linearly couples to the qubits at such bias condition. 

We suggest using an extra high-frequency qubit for 
the parametric AC-modulated coupling. The used circuit 
can be otherwise identical to the primary qubits, but the 
splitting should be larger so that the third qubit is always 
in the ground state. The functioning of this coupling can 
be viewed also in terms of the so-called quantum induc- 
tance. In fact the present scheme is very similar to us- 
ing a Cooper pair bojsi as a coupling element or using 
an RF-SQUID 15 both of which schemes were suggested 
as ordinary DC coupling. The use of a third qubit has 
many nice features. First of all, it should be easy to fab- 
ricate this kind of circuitry since the technology is com- 
patible. It is easy to ensure that the third qubit indeed 
has higher frequency than the other two. Secondly, the 
control can be all magnetic and frequency multiplexed 
(assuming different splittings) with a minimal number of 
microwave lines. Also, it turns out that the effective cou- 
pling at DC can be set to zero if desired. This however 
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may not be necessary because of the a x ® cr^-nature of 
the coupling. What is more important is that a kind 
of an optimal point exists for the effective coupling en- 
ergy. This may be crucial since at this point the effective 
coupling energy is insensitive to low- frequency flux noise, 
and therefore two-qubit oscillations can be expected to 
be long lasting. The time to generate a universal entan- 
gling gate similar to CNOT is on the order of 10-100 ns 
with experimentally realistic parameters. 

The paper is divided as follows. In Sec. [n] we first de- 
rive an approximate Hamiltonian starting from the three- 
qubit Hamiltonian. The meaning of the different terms 
will be discussed. Then in Sec. IIIII we discuss the imple- 
mentation of the coupling scheme in the context of the 
present system. In Sec. lIVI we present results of numerical 
simulations. There the validity of the effective Hamilto- 
nian is tested in simulation by comparing the dynamics 
of the truncated system with that of the full system. Sec- 
tion [V] is dedicated to discussion. 

II. SYSTEM AND THE HAMILTONIAN 

We consider a system of three inductively coupled flux 
qubits as illustrated in Fig. ^ The qubits consist of, say, 
three^ or four— ideally identical junctions except for one 
junction that has an area smaller by a factor a. The 
large junctions are characterized by a charging energy 
Eq = e 2 /2C where C is the junction capacitance and 
by a Josephson energy Ej = hlc/2e where Iq is the 
critical current of the junction. Ideally, the smaller junc- 
tion has a charging energy Ec/a and a Josephson en- 
ergy aEj. Typically, Ej/E c ~ 10 1 - 10 2 . For three- 
junction qubits a is around 0.75 and for four-junction 
qubits around 0.5. These choices ensure that a two-well 
potential forms and furthermore, a two-level (qubit) ap- 
proximation of the system is valid. We neglect the effect 
of loop inductances in the single qubit energies, which 
should be valid as long as Josephson inductance dom- 
inates. For a rigorous derivation of the coupled-qubit 
Hamiltonian, see Ref. |l9j. 

The three-qubit Hamiltonian reads in the persistent 
current basis 

H 3q = H 3 dc + H 3 mw (1) 
The DC part of the Hamiltonian reads 

1 3 1 
H 3DC = - 2 E + ~ o E ( 2 ) 

j=i k^i 

where tj = 2I p j(&j — $o/2), I p j is the maximum persis- 
tent current of the qubit j, <I>j is the flux threading the 
loop of the j th qubit and <I>o = h/2e is the flux quantum. 
The approximation is valid near the optimal point for the 
qubits, i.e., ^ ss <&o/2. The tunneling amplitudes Aj are 
defined by the Josephson energies and capacitances of the 
qubit junctions. These splittings depend in particular ex- 
ponentially on a. Moreover, the couplings are defined by 
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FIG. 1: Schematic illustration of three anti-ferromagnetically 
coupled flux qubits. 



the mutual inductances M kh i.e., J k i = Jik = M kl I vk I pl . 
The couplings can be either ferromagnetic (M k i > 0) or 
anti-ferromagnetic (M k i < 0) depending on the geome- 
try. The microwave controlled Hamiltonian simply reads 



H 3M w = - 2 E < M*) cr ^ ( 3 ) 
i=i 

Here 5tj(t) = 2I p j5<&j(t), and S®j(t) is the AC compo- 
nent of flux threading the loop of the qubit j. 

Let us assume that u>i,u>2 *C w 3 , where ujj = 

yj A 2 + e 2 . Applying the unitary transformation 

yields (H 3q = UH 3q W) 

1 2 1 

Hz q = - 2 E + e i°2) - 2^^ " Jl2fJ ^ 
j'=i 

+ - ^)"-) - \p'^< 

(5) 

We denote the identity operation on qubit j by Ij . Here 
we have included the temporal dependence of the bias of 
the third qubit in the transformation, i.e., we are trans- 
forming to the adiabatic basis of the third qubit. That 
is e 3 (t) = e 3 + 5e 3 (t) and uj 3 {t) = ^fe 3 (t) 2 + A 2 . We will 
use an adiabatic approximation for the third qubit. 

By looking at Eq. (JHJ) we see that in the lowest order 
adiabatic approximation for the third qubit, i.e., when 
o\ — > (trf) = 1 and a\ — > (cr^) = 0, the flux biases of the 
qubits 1 and 2 will be shifted. The c^-component will 
have no effect on the dynamics of the two qubits to low- 
est order. But going to higher order in J13 and J23 will 
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yield an effective coupling term. In order to eliminate and 

the third high-frequency qubit we use a trick known as 

the Schrieffer- Wolff transformation (see e.g. Ref. pcj)- 

We look for a transformation exp(— S) such that the an- 

tihermitian operator S = — S' is first order in J^i such 

that it eliminates the cr^-component in first order. That P 

is, we require 



[S,H ) = H 1 +O(Ji l ) 



(6) 



J13A1A3W3 



J23A2 A3CJ3 2 

(u,3to 2 -Ai)Mt) 2 -ei)-Al e l <T »- 



(14) 



where 
#o = 
and 



1 2 1 

-22 (Ajcr^. + e,o^) - -w 3 (t)<r! - Ji 2 cr]cr 2 (7) Now the effective Hamiltoman reads 



#1 = ( + J 23 CT 2 ) 



A3 3 



(8) 



'ws(t) *' 

If Eq. © is satisfied we get up to second order 
exp(-S)(H Q +H 1 )exp(S) 

= H + H 1 + [H Q +Hi,S] + ±[[H ,S},S] + 0(4) 

= H + (Hi - [H , S]) +~[H 1 ,S\ + 0(4) (9) 

such that the term Hi is eliminated to first order since 
(Hi — [Ho,S]) is second order. The total Hamiltonian 
has also the terms 



cxp(-S)(ff + H X +H 2 + H 3 ) exp(S) 
*Ho + ~[Hi,S}+H 2 + H 3 . 



(15) 



The second order terms [H 2 , S] and (Hi — [Ho, S]) were 
dropped because they are off-diagonal in the qubit 3 op- 
erators. Since the single-qubit MW-term H 3 is anyway 
weak, we will not consider its transformation. We also 
neglect dS/dt since it vanishes when the adiabatic ap- 
proximation is made. We need to calculate 



and 



H 2 = -(Ji 3 al + J 23 a 



2\^W ^3 

W 7 * 



1 2 



3 = 1 



(10) 



(11) 



2 x u 3 (t) 



5 (Jl3 a l + ^23f 2 ), 



(16) 



It turns out that Eq. JJJ is solved by choosing 



S = 



(12) 



where {• , •} stands for the anticommutator. We get 



where a and j3 contain only operators operating on the 
qubits 1 and 2. A straightforward calculation gives up to 
first order 

J13A3 1 



+ 



u 3 (ty - a\ - 

J2 3 A 3 



J13A3A161 1 

(o; 3 (t) 2 -Af)(a;3(t) 2 -e?)-A?e?^ 
J23A3A262 2 



(13) 



2 WO 



(Jl 3 oi + J23Cr 2 ),^] = 

^AiA? , 



u> 3 (i) 2 M*) 2 - A? - e 2 ) * 



+ 



J 23 A 2 A§ 



W3 W 2 (^3(i) 2 -A| 



el)^ 



(17) 



and (up to constant when the approximation w 1 is 
made) 
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2 i u 3 (t) 



(J 13 crl + J 23 a 2 ),a}a 3 z 



u 3 {tf 



J23 J13A3 
J23J13AI / A l£l 



4 



A 2 
^2 



A 2 - e 2 
1 e i 

A 2 e 2 



A 2 
^2 



(18) 



We may safely neglect the cr^cx 2 and a^cr 2 terms if the 
third qubit frequency u 3 (t) is significantly higher than ej 
and Aj for j = 1,2. The effective Hamiltonian for two 
qubits can thus be compactly rewritten as 

1 2 1 2 

H *« = - 2 £ + g ^0 -JuMtt- 2 E m*H 



3=1 



3=1 



with 



e 1 = e 1 



2e 3 (t)J j3 
W3(t) ' 



2(J J3 A 3 ) 2 A J 



w 3 (f) 2 Mt) 2 



A 3 2 



(19) 
(20) 

(21) 



and 



= J: 



12 



J23 J13A 2 
w 3 (t) 3 W*) 



wa(*) 2 



w 3 (t) 2 - ef 



A 2 
^2 



w 3 (t) 2 



A 2 



(22) 



Let us discuss the physical meaning of the different terms 
in the approximate Hamiltonian in Eq. 119|) . Recalling 
that Jj 3 = Mj 3 I p jI p3 and since the ground state expecta- 
tion value of current for the third qubit is in the absence 
of the coupling 



<0 3 |/ 3 |0 3 ) = -i^-^3 : 
2 a® uj 3 



(23) 



we see immediately that the correction in Eq. ill' 01) 
2e 3 (t)Jj 3 /u 3 (t) = —2Mj 3 I p j(0 3 \I 3 \0 3 ) can be interpreted 
as a shift of the qubit j bias due to the circulating cur- 
rent of the third qubit. From now on we consider the 
modified "sweet spot" ej = 0. Note that this cannot 
be perfectly achived at DC because the time average of 
?j = depends on the modulation amplitude of the third 
qubit microwave. However, the deviation from the ideal 
case is only second order in <5e 3 (t) / io 3 , where oj 3 is the DC 
value of the third qubit frequency and 5e 3 (t) is the mi- 
crowave modulation of the third qubit bias. Other than 
that, we neglect the time dependence of ij since it can 
be absorbed in the microwave Hamiltonian and results 
in a weak crosscoupling. Since we are considering a res- 
onant control scheme such omissions are well justified. 



The interpretation of the renormalization of Aj is more 
complex. 

The effective coupling term in Eq. (|22|) can most easily 
be understood as coupling through the quantum induc- 
tance of the third qubit. By quantum inductance we 
mean the inverse of the coefficient by which the circulat- 
ing current of the auxiliary qubit responds to a change 
in its bias flux, i.e. 



1 



L Q (t) 



d(0 3 \h\0 3 ) 



2/ 2 3 Aj 
^ 3 W 3 ' 



(24) 



In the limit of large uj 3 we therefore obtain 



Ji 2 (*) » J 12 + 



2 J2 3 </i 3 Ao 



Mis 



M 23 Mi; 

L Q {t) 



p2- 



(25) 



From this form it is clear that the effective mutual in- 
ductance between the qubits is affected by the quantum 
inductance of the auxiliary qubit; qubit 3 generates a 
shielding current in response to changes in the flux of 
qubit 1 (2) which in turn changes the flux through qubit 
2 (1). In order to cancel the effective DC coupling the 
following conditions must be satisfied since Lq < 0: If 
Mi 2 < 0, then Mi 3 and M 23 must have the same sign 
but if Mi 2 > then Mi 3 and M 23 must have a differ- 
ent sign. Naturally the magnitudes need to be suitable, 
too. Figure |2 illustrates the coupling term as a function 
of flux with realistic experimental parameters. Note that 
two points exist where go = 0. 

To conclude the present Section we rewrite the Hamil- 
tonian conveniently at the optimal point £j =0, via ro- 
tating it by exp(i7r/4(CT^ + er 2 )), as 



rropt 
n 1q 



1 2 

- £ (h 3 a\ - Se 3 (t)ai) - Ji 2 (t)<r> 2 



3 = 1 



(26) 



Because at the optimal point observed single qubit co- 
herence times are much superior to those measured else- 
where, we will focus our attention there. 



III. IMPLEMENTATION OF THE 
PARAMETRIC COUPLING SCHEME 

As suggested in Ref. , the form of the Hamiltonian 
in Eq. l|26|l is ideal for the realization of a coupling scheme 
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of the microwave phase. All the control is carried out by 
a combination of phase and amplitude modulation. We 
have set the phase of u>± to zero for simplicity. 

The Hamiltonian in Eq. (|3L)[1 is clearly universal. All 
the single qubit operations can be obtained easily e.g. 
by rotating three times around two different axis. From 
the point of view of two-qubit operations Eq. (|3U|) is ex- 
tremely nice because it can in principle be used to realize 
directly the so-called B-gate2i, which is known to be the 
best all-purpose two-qubit gate there is. Namely two ap- 
plications of the B-gate suffice to realize any two-qubit 
gate. Essentially, the B-gate is up to some local one-qubit 
operations given by 



0.45 0.46 0.47 0.48 0.49 0.5 0.51 0.52 0.53 0.54 0.55 
3 



Ub ~ exp 



(7 (7 



2 a v a v 



(31) 



To realize this with a MW pulse of the length At one just 



FIG. 2: The dependence of the effective coupling term on sets Qj = 0, g+At/h = 1/8 and g_At/h = 3/8. 



flux when I pj = 0.4/xA, A 3 /h = 10 GHz e 2 /h = ei/h = 1.72 
GHz, A 1 /h = 3.1 GHz, A 2 /h = 4.3 GHz, , Mia = M 23 = -6 
pH and M12 = — 1 pH. (For these choices ei = 62 = if 
e 3 /h = 7.375 GHz.) 



in which the coupling constant J\2(t) is modulated sinu- 
soidally at the angular frequency u>± — (A 2 ± A{)/K. 
The essence of the scheme is seen easily by considering a 
general modulation of the form 

Juit) = g + g+(t)cos(uj + t) + cos(w_t). (27) 

Here 



dJ 



g± 



12 



de 3 



Sc. 



3± 



(28) 



where 8e^± is the amplitude of the modulation of e$(t) at 
either the sum or the difference frequency and go is the 
DC component of the coupling Jn{t) of Eq. (|22|) . For 
single-qubit operations we use Rabi oscillations driven 
by a resonant microwave 



Sej(t) = 2ttj(t) cos(Ajt/h + (j>j(t)). 



(29) 



In our setup all the temporal dependence of the Hamil- 
tonian is assumed to arise from the time-dependent flux. 
The rotating wave approximation, which is also valid if 
crosscouplings are taken into account, results in a rotat- 
ing frame Hamiltonian of the form 



HI 



1 2 

^ = 5 Efli(t) (cob 



9+(t) 



9-(t) 



- sin 



^„ 2 \ 



(30) 



This can be taken as our logical Hamiltonian. The ap- 
proximation is valid when go,g± < w± and fL < A,-. 
We must include the possibility of phase shift to get ar- 
bitrary one qubit gates, hence the temporal dependence 



If the modulation of the coupling is not possible at 
both the sum and the difference frequencies, say, due to 
a too small difference in lo + and W3, we may still do a 
universal quantum gate when g + — 0. The easiest way 
to achieve this is to do a gate essentially equal (i.e. locally 
equivalent) to the so-called double-CNOT, or DCNOT. 
The nontrivial two qubit part of this gate can be written 
e.g. as 



(7T 
*4 + <r v °y) 



(32) 



Choosing flj = 0, g+ = and g-At/h = 1/2 realizes the 
gate. Three applications of J7dcnot with some one-qubit 
gates suffice to realize any two-qubit gate2£. The gate is 
thus just as good as the better known CNOT. 

In order to seriously evaluate the viability of the 
present scheme, one needs to consider the matter of de- 
coherence. We are operating in the optimal point of the 
individual qubits so that there is more hope in this re- 
spect than in some non-optimal point scenario. However, 
in the presence of coupling, the situation is more compli- 
cated. As has been previously showni^, if go <§; w_ , the 
coupling circuit does not add seriously to decoherence, 
meaning the single-qubit transition frequencies are sta- 
ble. Especially if go — 0, the individual qubits are truly 
first order insensitive to fluctuations in all the fluxes Qj. 
If go is however too large compared to uj- then the ben- 
efits of the optimal point may unfortunately be lost. 

It is also important to consider the phase coherence of 
the |01) to |10) (or |00) to |11)) oscillation. It turns out 
that our scheme is nice in this sense. Namely, we would 
like the frequency of oscillation to be stable against low- 
frequency noise, i.e. dg±/d$3 — 0. There indeed exist 
points symmetrically around $3 = $o/2 at which this is 
satisfied, i.e., approximately at £3 = ±As/2. Figure |3| 
illustrates the form of g± with the same parameters as 
in Fig. |2 However, there is no reason for these spe- 
cial points to be ones in which go = 0. For perfect first 
order insensitivity to 1/f flux noise, it is necessary to 
design the circuitry so that the conditions go = and 
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FIG. 3: The effective AC coupling constant g± as a function 
of DC flux when the flux amplitude through the third qubit 
has the realistic value <5$3 — 2 x 10 -4 $o- 



dg±/d&3 = are fulfilled simultaneously. Currently the 
technology may not be there yet though, since in partic- 
ular A3 is hard to control in fabrication. It should be 
however possible to satisfy go <C o-i- well enough and set 
still dg±/d&3 = 0. One possible solution is to fabricate 
qubits with comparably high Aj . 

Let us compare the proposed scheme with the ear- 
lier scheme of coupling through a DC-SQUID. In our 
scheme the added improvement is the existence of the 
additional stability in the form of dg±/d&3 = 0. Also, 
the bias conditions are more favorable for the individ- 
ual qubits. Although the effect of the DC-SQUID on 
the single qubit terms analogous to Eqs. H20I21JI was not 
included in the consideration of coupling through a DC- 
SQUID in Refs- HHj , the SQUID will couple directly to 
the individual qubits. At high bias current and close to 
$o/2 the bias current noise always couples in first order 
to <j x . This is known to especially cause a rapid decrease 
of the relaxation time T-^ and thus also of the dephas- 
ing time. Because of the expected significant decline in 
Ti we argue that the use of a third high-frequency qubit 
may be a more viable option for realizing the promising 
parametric coupling scheme. 



IV. SIMULATION 

In order to confirm the validity of the approximations 
made, we present some simulation results. We can con- 
veniently simulate the dynamics of the system both with 
the full three-qubit Hamiltonian and with the effective 
Hamiltonian in Eq. (|26|l . The form of the full Hamilto- 
nian that we use is the Hamiltonian in Eq. Q rotated 
by the time-independent version of Eq. J3J so that the 
Hamiltonian is exact but in a more convenient basis. Also 



the qubits 1 and 2 are rotated. That is, 
1 2 1 

- Ji2el&l + (Ji3<rl + J23&1) [ — o\ -al ) 

\uj 3 uj 3 ) 

+ jE-MtK ~ \ {-^ ~ M*)- (33) 

Z ' Z \U>3 W3 J 

No rotating wave approximation has been used for the 
simulation, but rather we solve the full Schrodinger equa- 
tion. The results are obtained using Floquet states^ 
since the problem is time-periodic for the duration of the 
two-qubit pulse. Below, the initial state of the qubit 3 is 
always assumed to be IO3) in the basis of the above Hamil- 
tonian. The reduced density matrix p\^. is obtained by 
tracing over the third qubit. We concentrate here on 
the use of the difference frequency since this should be 
experimentally easier to achieve. 

As an illustrative example, let us consider the initial 
state I L1O2) • Figure 0] illustrates the probability of dif- 
ferent outcomes in the single-qubit eigenbasis vs. time 
with the example parameters of Fig. |2I 24 In this example 
the DC component of the coupling is go rs 0. In this 
simulation we assumed that the microwave couples only 
to qubit 3. However, even if the microwave crosscouples 
to all the qubit loops, it is still possible to carry out the 
operations with high fidelity. We see that after a time 
h/{2gJ) the probabilities of the states II1O2) and |0i I2) 
are flipped. This corresponds to a gate locally invari- 
ant with the DCNOT, i.e. the operation in Eq. 
The fact that the full three-qubit simulation shows re- 
duced amplitude after tracing over the third qubit is 
attributable to the fact that the third qubit entangles 
slightly with the primary qubits due to its finite fre- 
quency. At best the fidelity F — \J (V'|/>i.2|'0) between 
the reduced density matrix and the target pure state 
\tp) = exp(i7r/4(cr>2 +cr i CT 2^| li0 ^ ig about gg% _ 

The above example is not perfect but indicates clearly 
that the scheme should be comparably easy to demon- 
strate in an experiment. To see how increasing A3 af- 
fects the situation, we consider the case A3 = 15 GHz in 
Fig. [SJ It can be seen clearly that while the gate time 
stays reasonable (about 60 ns), the amplitude of the os- 
cillation and thus the fidelity is improved above 99%. 

So far we have given examples of operation at go w 0. 
It is however possible to work with finite go- Partic- 
ularly, we would like to operate at dg±/d$>3 = 0. To 
demonstrate, we simulate the first example again at this 
point. The results are shown in Fig. |SJ We now work 
at the frequency <D_ = oj- + 2<?q / cj_ . If this correction^ 
is not taken into account, the amplitude of oscillation 
is roughly halved but with the corrected operation fre- 
quency, the best fidelity is again about 99%. This is 
extremely promising since at this point the decoherence 
is expected to be reduced. 
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FIG. 4: (Color online) Example of the coherent oscillations 
induced by a microwave at the frequency (A2 — Ai )/h. Above 
the two-qubit approximation was used whereas the lower plot 
is the result of the full three-qubit simulation. The black 
line starting from probability 1 at t = corresponds to 
1 1 1 02 ) where as the blue line starting from probability at 
t — corresponds to IO1I2). The green and red lines cor- 
respond to the states 1 1 1 12 ) and 0i O2 ) , respectively. Here 
Ai ~ 27.5 ns results in gate equivalent to DCNOT. We use 
e-i/h = 7.375GHz. Here the microwave amplitude was #£3 = 
497 MHz which corresponds to Fig. [3] With these choices 
go ?s but dg±/d$ 3 0. 



Since we are interested in creating a unitary gate, it is 
not sufficient to consider one initial state alone. One pos- 
sible measure is to calculate the fidelities for a complete 
set of initial states. It however turns out that the fideli- 
ties with the other logical qubit states for the present 
scheme are equally good as for the state | I1O2) - This 
means that we can indeed carry out the DCNOT with 
high fidelity. 

We conclude that the present scheme is an experimen- 
tally realistic scheme for realizing two qubit gates at the 
optimal points of the qubits. We can also conclude that 
the derived two-qubit Hamiltonian is in good agreement 
with the full three-qubit Hamiltonian. The desired be- 
havior of the three-qubit system is only observed when 
the renormalization of e and A are taken into account. In 
the case of a finite go the operation frequency is further 
renormalized as can be seen in the example above. The 
predicted gate time also agrees reasonably well with the 
full dynamics. 

V. DISCUSSION 

We have shown how to realize a tunable coupling 
scheme for optimally biased flux qubits using an extra 
qubit as a coupling element. Although the mutual in- 
ductances are fixed, the primary qubits can be effectively 



FIG. 5: (Color online) Same as in Fig. 0] but with A 3 = 15 
GHz and e 3 /h = 5.73GHz. The fidelity is at best above 99%. 
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FIG. 6: (Color online) Evolution of the probabilities at the 
"sweet spot" of the coupling dg±/d§z — 0. The parameters 
are like in Fig. 0] otherwise, but t3/h = 4.7 GHz, ei = €2 = 
1.23 GHz and g /h = 110 MHz. Thus again £1 = £2 = 0. The 
fidelity is at best above 99% and the gate time is about 26 ns. 



coupled and decoupled at will using microwaves at a suit- 
able frequency. As a first test of the scheme, we suggest 
the coupling to be demonstrated with just one measure- 
ment SQUID with all the three qubits sitting inside the 
loop. The seeming visibility problem between |0i I2) and 
1 1 1O2) in this kind of a measurement can be overcome 
easily. To realize coherent two-qubit oscillations of the 
kind presented in the previous Section and to see a sig- 
nature of the tunable coupling, one can convert the oscil- 
lation in the subspace spanned by IO1I2) and | I1O2) into 
an oscillation of IO1O2) and 1 1 1 12) by just one 7r-pulse. 
Scanning the MW frequency around w_ should give ad- 
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equate proof of the tunability. This experiment could be 
carried out even with a continuous wave at w_ and two 
chopped 7r-pulses, one to prepare either |0i I2) or | I1O2) 
and one to overcome the visibility problem. Scanning the 
interval between the pulses should result in an oscillation 
between IO1O2) and 1 1 1I2) - 

For more sophisticated control of the qubits, all the 
microwaves need to be chopped and attention needs to 
be paid to phases. It is of course important that the sig- 
nals are phase-coherent. The choices made in this paper 
should be realizable in an experiment by properly ad- 
justing the phase of the microwave source addressing the 
qubits. Once this is done the exact timing of the coupling 
pulse is not important as time does not appear explicitly 
in the rotating frame Hamiltonian. In the case of the 
FLICFORQ scheme timing is more crucial since the ro- 
tating frame Hamiltonian still has an explicit time depen- 
dence. The phases in Eq. (|3L)[) could of course be chosen 
differently from what we have done. Like in FLICFORQ, 
it should be possible to utilize frequency multiplexing in 
the present coupling scheme. 

Let us comment on the scalability of the present 



scheme. A generalized version of the presented transfor- 
mation may be used to derive an effective Hamiltonian 
for a larger multiqubit register without changing the form 
of the result. For instance in the case of a linear chain 
of qubits with every second qubit acting as a coupling 
element the bias <Lj of any computational qubit will have 
a contribution from both of its neighboring coupling el- 
ements and the same holds for Aj. The coupling be- 
tween neighboring computational qubits will be exactly 
the same as in the text up to second order. We can there- 
fore conclude that the scheme is potentially scalable to 
multiple qubits since figuring out the controls required for 
a more complicated multiqubit unitary gate is a tractable 
task. 
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